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Abstract 

We discuss a general class of nonlinear mean-field Fokker-Planck equations [P.H. Chavanis, Phys. Rev. E, 68, 
036108 (2003)] and show their applications in different domains of physics, astrophysics and biology. These 
equations are associated with generalized entropic functionals and non-Boltzmannian distributions (Fermi-Dirac, 
Bose-Einstein, Tsallis,...). They furthermore involve an arbitrary binary potential of interaction. We emphasize 
analogies between different topics (two-dimensional turbulence, self-gravitating systems, Debye-Hiickel theory of 
electrolytes, porous media, chemotaxis of bacterial populations, Bose-Einstein condensation, BMF model, Cahn- 
Hilliard equations,...) which were previously disconnected. All these examples (and probably many others) are 
particular cases of this general class of nonlinear mean-field Fokker-Planck equations. 



1. Introduction 

Recently, several researchers have questioned the "universality" of the Boltzmann distribution in physics. This 
problem goes back to Einstein himself who did not accept Boltzmann's principle S = k\nW on a. general scope 
because he argued that the statistics of a system {W) should follow from its dynamics and cannot have a universal 
expression [1]. In 1988, Tsallis introduced a more general form of entropic functional in an attempt to describe 
complex systems [2]. This was the starting point for several generalizations of thermodynamics, statistical me- 
chanics and kinetic theories. A lot of experimental and numerical studies have then shown that complex media 
exhibit non-standard (non-Boltzmannian) distributions and that Tsallis q-distributions can fit a wide diversity 
of results when the Boltzmann distribution fails. However, there also exists situations that are described neither 
by the Tsallis nor by the Boltzmann distribution. The important point is to explain why a given system exhibits 
non-Boltzmannian distributions and find the reason for that. We cannot simply invoke Tsallis entropy each time 
the system is "complex" . 

In a recent series of papers [3], we have argued that the Tsallis distributions and the associated entropic 
functional are just particular (yet important) types of non-Boltzmannian distributions and entropies, and that 
there is no fundamental reason why they should play a prominent, or universal, role in physics. Other distributions 
can emerge as well in complex systems for different reasons. On general grounds, we have suggested that non- 
Boltzmannian distributions arise when "hidden constraints" are in action and prevent the a priori accessible 
microstates to be equiprobable. These constraints can be of different nature: (i) There can be exclusion or inclusion 
constraints on the probabilities of transition of the system from one state to the other. A fundamental example is 
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the Pauli exclusion principle in quantum mechanics which prevents two fermions to occupy the same site in phase 
space. This leads to the Fermi-Dirac statistics, instead of the Boltzmann statistics, placing a bound /(r, v, t) < r?o 
on the distribution function. Similarly, one can account for excluded volume constraints in classical systems of 
particles. In simplest models, this loads to the Formi-Dirac statistics in physical space placing a bound p{r, t) < ao 
on the spatial density. Alternatively, we can account for inclusion constraints (like for bosons) by favoring the 
transition to a site that is already occupied. This leads to the Bose-Einstein statistics. Nonlinear Fokkcr-Planck 
equations taking into account these inclusion or exclusion constraints have been discussed by Kaniadakis [4] . (ii) We 
can imagine situations where some microscopic constraints prevent the system from sampling the energetically 
accessible phase space uniformly. These microscopic constraints will modify the equilibrium distribution with 
respect to the expected Boltzmann distribution. One example is the case of porous media where the system has a 
complicated phase space structure (fractal or multi-fractal) leading to anomalous diffusion. The TsaJlis entropy [2] 
is adapted to this situation. The dynamical evolution of the system can be described by a nonlinear Fokker-Plamck 
equation introduced by Plastino & Plastino [5]. (iii) In biology, the diffusion coefficient and the mobility of the 
particles arc often assumed to depend on the local concentration so as to account for microscopic constraints 
that cannot be easily modeled. This is the CEise for the Keller-Segel model of chemotaxis [6] which can lead to 
non-standard distributions at equilibrium, (iv) Finally, coUisionless stellar systems and inviscid 2D turbulent flows 
described by the Vlasov-Poisson system and by the 2D Euler-Poisson system spontaneously form meta-equilibrium 
states (galaxies and vortices) that are usually described by non-Boltzmannian distributions [7]. This is due, on 
the one hand, to the existence of an infinite family of Casimir invariants for the fine-grained dynamics which lead 
to non-standard distributions on the coarse-grained scale. These Casimirs play the role of "hidden constraints" 
because they are not directly accessible from the coarse-grained flow where the observations are made. Generalized 
Fokker-Planck equations have been introduced by Robert & Sommcria [8] and Chavanis et al. [9] to take into 
account these constraints. On the other hand, coUisionless mixing is incomplete in general and other forms of 
distributions can emerge that arc different from those predicted by the statistical mechanics of violent relaxation. 
These quasi-stationary states (QSS) are nonlinearly dynamically stable stationary solutions of the Vlasov or the 
2D Euler equations selected by the dynamics. In this context, the Tsallis distributions are particular stationary 
solutions of the Vlasov equation known as stellar polytropcs in astrophysics. They can sometimes arise as a result 
of incomplete violent relaxation, like in 2D turbulence, but this is not generic [10]. For coUisionless stellar systems 
and 2D vortices described by the Vlasov and the 2D Euler equations, non-standard distributions are explained by 
microscopic constraints (the Casimirs) and incomplete violent relaxation (non-ergodicity) [11]. 

In this paper, we discuss a general class of nonlinear mean-field Fokker-Planck equations associated with gen- 
eralized ontropic functionals (loading to non-standard distributions) and an arbitrary form of binary potential of 
interaction [3]. These equations are consistent with an effective generalized thermodynamical formalism (E.T.F.). 
The study of these equations is interesting in its own right due to their rich mathematical properties but we provide 
here explicit examples (2D turbulence, self-gravitating systems, porous media, bacterial colonies, Bose-Einstein 
condensation...) where these nonlinear mean-field Fokker-Planck equations can have physical applications. 



2. Nonlinear mean-field Fokker-Planck equations 

We consider a generalized class of nonlinear mean- field Kramers equations of the form [3]: 



where / = f{r,v,t) is the distribution function, C is a convex function (i.e. C" > 0), /3 = 1/T is the inverse 
temperature and <I>(r, t) is a mean-field potential. The friction parameter satisfies an Einstein relation ^ = D(5. 
Since the temperature T is fixed, the above equation describes a canonical situation. We assume that the potential 
is related to the density by a general relation of the form 




(1) 




(2) 



where w is a binary potential of interaction and p = J fdv is the spatial density. We introduce the functional 
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This functional can be interpreted as a free energy F = E — TS, where 5* is a "generalized entropy" and E — K+W 
is the energy including a kinetic term and a potential term. When S[f] = — J / In fdvdv is the Boltzmann entropy, 
Eq. (1) reduces to the ordinary Kramers equation. More generally, it is straightforward to check that F[f] plays 
the role of a Lyapunov functional satisfying _F < 0. This is the equivalent, in the canonical ensemble, of the 
H-theorem. The stationary solutions /^^(r, v) are determined by the integro-differential equation 



C"(/e,) = 



+ *(r) 



(4) 



where 3> depends on feq by Eq. (2). This equation results from the individual cancelation of the advective term 
(L.H.S.) and collision term (R.H.S.) in Eq. (1) [3]. Since C is convex, Eq. (4) can be inversed so that feq = feq{() 



is a decreasing function of the individual energy 



/2 + $(r). The linearly dynamically stable stationary 



solutions of Eq. (1) arc rmmrna of the free energy F[f] at fixed mass [3]. 

In the strong friction limit ^ +oo, or for large times t ^ the distribution function /(r, v,t) is close to 
the isotropic distribution function determined by the relation 



C'if) = -P 



+ A(r,t) 



+ o(r')- 



(5) 



The pressure p = ^ffv^dv = p[\{r,t)] and the density p = f fdv = p[X(r,t)] are related to each other by a 
barotropic equation of state p = p{p) which is entirely specified by the function C(/). This equation of state 
is obtained by eliminating A(r, t) between the foregoing relations. Taking the hydrodynamic moments of the 
generalized Kramers equation (1) and closing the hierarchy by considering the limit g — > +oo, it is shown in [3] 
that the time evolution of the density p(r, t) is governed by the generalized Smoluchowski equation 



dp 
dt 



= V- 



:(Vp + pV$) 



(6) 



A formal derivation of the generalized Smoluchowski equation (6) from the generalized Kramers equation (1) can 
be realized by performing a Chapman-Enskog expansion in powers of [12]. The generalized Smoluchowski 
equation monotonically decreases the Lyapunov functional 



p 

F[p\ = j pj ^dp'dv+\ I p$dr. 



(7) 



which is the simplified form of free energy (3) obtained by using the isotropic distribution function (5) to express 
F[f] as a functional of p [3,12]. A stationary solution satisfies the condition of hydrostatic equilibrium 



Vp + pV$ = 0. 



(8) 



Taking the hydrodynamic moments of the nonlinear mean-field Kramers equation (1) and closing the hierarchy by 
a condition of local thermodynamic equilibrium (in the canonical ensemble) one gets a system of hydrodynamical 
equations 



dp 
dt 



+ V • (pu) = 0, 



^ + (u ■ V)u = -ivp - V3- - Cu, 
at p 



(9) 
(10) 



that we called the damped barotropic Euler equations [3]. In the strong friction limit ^ — > +oo, one can neglect 
the inertial term in Eq. (10) to obtain pu = — |(Vp-|- pV$) + 0{^~^). Substituting this relation in the continuity 
equation (9), we recover the generalized Smoluchowski equation (6). The generalized Smoluchowski equation can 
also be written 



dp 
dt 



pC"{p)Vp + 0pV^ 



(11) 
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In the general formulation, D can depend on r and t, so on p(r, t) [3]. The generalized Smoluchowski equation 
(11) monotonically decreases the generalized free energy 




(12) 



The stationary solutions of Eq. (11) are solutions of the integro-differential equation 



C'(Pe,) 



/3$(r) - a, 



(13) 



where $ depends on p^q by Eq. (2). Since C is convex, Eq. (13) can be inversed so that peg — peg ($) is a decreasing 
function of the potential (assuming (3 > 0). The linearly dynamically stable stationary solutions of Eq. (11) are 
minima of the free energy F[p] at fixed mass [3]. 

These generalized Kramers and Smoluchowski equations can be obtained in different ways: (i) from the linear 
thermodynamics of Onsager [13,3] (ii) from a variational principle called Maximum Entropy Production Principle 
(MEPP) adapted to the canonical description [3] (iii) from stochastic Langevin equations where the diffusion 
coefficient (noise) and the friction or mobility coefficients explicitly depend on the density of particles [14,3] (iv) 
from the Master equation, by allowing the transition probabilities to depend on the population in the initial and 
arrival states [4] . These generalized kinetic equations can be justified as effective equations modeling heuristically 
microscopic constraints that arc not directly accessible to the observer. It is clear that the mathematical study of 
these equations is of considerable interest. A systematic study of the generalized mean-field Smoluchowski equation 
(6) has been undertaken in [15,16,17,18,19,20] for different equations of state p(p), or generalized entropies S[p], 
and different potentials of interaction M(r — r') in different dimensions of space d. Note that when the potential 
of interaction is short-range, the mean-field potential is given by $ ~ —p — R^Ap and, when substituted in Eq. 
(11), wc get a generalized form of the Cahn-Hilliard equation (see [3,12]). We now give explicit examples where 
these nonlinear mean-field Fokker-Planck equations have physical applications. 



3. Two-dimensional turbulence 

Two-dimensional flows with high Reynolds numbers are described by the 2D Euler equations 



In geophysical fluid dynamics, one also considers the QuEisi-Geostrophic (QG) equations where the vorticity to is 
replaced by the potential vorticity q related to the streamfunction ip hy q = —Atp + where R is the Rossby 
radius. Starting from a generically unstable initial condition uo(r), the 2D Euler equations are known to develop a 
complicated mixing process which ultimately leads to the emergence of a large-scale coherent structure, typically 
a jet or a vortex. Jovian atmosphere shows a wide diversity of structures (Jupiter's great red spot, white ovals, 
brown barges, ...). One question of fundamental interest is to predict the structure and the stability of these 
"equilibrium" states depending on the initial conditions. To that purpose, Robert & Sommeria [21] have proposed 
a statistical mechanics of the 2D Euler equation. The idea is to replace the deterministic evolution of the flow 
tj(r,t) by a probabilistic description where p{r,a,t) gives the density probability of finding the vorticity level 
w = 0" in r at time t. The observed (coarse-grained) vorticity field is then expressed as uj{r, t) = J pada. 

Consider first the case where the initial vorticity field consists of an ensemble of patches with vorticity u) = ao 
surrounded by irrotational flow a; = 0. These patches will mix in a complicated way but their vorticity and their 
total area 70 will remain constant. In this two-levels approximation, the conservation of the area is equivalent to 
the conservation of the circulation P = J Zodr = 7ocro. The energy E = ^ J ujipdr is also conserved. Using a result 
of large deviations, it can be shown [21] that the optimal probability p(r) of finding the level ctq in r at equilibrium 
is obtained by maximizing the mixing entropy 



at fixed 70 = / pdr (or F) and E. This expression of entropy can be obtained from a combinatorial analysis taking 
into account the "incomprcssibility" of the vorticity field. Using the method of Lagrange multipliers, the optimal 
probability p*{r), and consequently the coarse-grained vorticity ZJ = p,ao, is given by 




u = — z X Vtp, 



uj = -Atp. 



(14) 




(15) 
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^i^)^-. ,,, , (A>0). (16) 

Note that aJ < ao, as the coarse-grained vorticity locally averages over the levels uj = and lj — ao. This constraint 
is similar to the Pauli exclusion principle in quantum mechanics and this is why (15) and (16) look similar to the 
Fermi-Dirac entropy and the Fermi-Dirac distribution. Non-standard (non-Boltzmannian) distributions arise here 
because the patches (levels) cannot overlap (or cannot be compressed). 

Robert & Sommeria [8] have proposed a dynamical model to describe the relaxation of the system towards 
statistical equilibrium. By using a Maximum Entropy Production Principle (MEPP), they obtain in the two-levels 
case an equation of the form 

^-hu-Va7= V-[D(Va7 + /3(t)a;(cro-u;)V^/>)], uJ = -Aip, (17) 

^(^^--/^fe^^^KV^' ^('•.*)«--^(-o-^). (18) 

This relaxation equation enters in the class of generalized mean-field Fokker-Planck equations of Sec. 2. It is 
associated with the Fermi-Dirac entropy in physical space S\uj] — — J {(ZZJ/cro) ln(ZZ;/(Jo)-|-(l— U/ao) ln(l— ZZ;/ao)}dr 
and with a potential it(|r — r'|) = — In |r — r'|. Note, however, that in the present case the inverse temperature 
f3{t) evolves with time according to Eq. (18) in order to conserve the energy. Therefore, Eqs. (17)-(18) increase 
the entropy S at fixed E and P. This corresponds to a microcanonical description while the equations presented 
in Sec. 2 correspond to a canonical description [3]. 

In the general case, the optimal vorticity distribution p(r, a) maximizes the mixing entropy 

S[p] = - J plnp drda, (19) 

at fixed E, F and F„>i — J uj"dr — J pa"dadr. The last constraints are equivalent to the conservation of the 
Casimirs or of the area 7(0-) = J pdr of each vorticity level. We must also account for the local normalization 
condition J pda — 1. This yields the Gibbs state 

+00 

Pi^^'')-^xi'^)e-^^''^"^", Z= j x(a)e-(^'*+'^)-da, (20) 

— 00 

where xi'^) ~ exp(— X]n>i Q^"""") encapsulates the Lagrange multipliers a„>i accounting for the conservation of 
the fragile constraints r„>i while /3 and a account for the conservation of the robust constraints E and F {Z is 
the normalization constant determined by J pda = 1 The coarse-grained vorticity uj = J pada is now given by 

where F{^) — — (lnx)'("I') and x(*&) = x{'^)^~'^^d,a. In this approach, xi*^) is determined a posteriori by the 
initial conditions through the conservation of the Casimirs. Indeed, the a„ are Lagrange multipliers which must 
be related to the constraints r„. The general relaxation equations have the form 



dt 



+ u-Vp = V ■ [D{Vp + P{t)p{a - iU)VV')], D (XL02^ J p{a-l 



(22) 



where the time evolution of P{t) is determined by the energy constraint [8]. These equations for p(r, a, i) conserve 
all the Casimirs (in addition to the circulation and the energy) and increase the mixing entropy (19) until the 
maximum entropy state (20) is reached. Note also that the diffusion coefficient is not constant but depends on 
the local fluctuations of vorticity. This can block the relaxation in a sub-region of the flow (maximum entropy 
bubble) and account for incomplete relaxation [8,9]. 



^ E and T are called robust constraints because they can be expressed in terms of the coarse-grained field. By contrast, the 
rTi>i arc called fragile constraints because they are not conserved by the coarse-grained vorticity as uj'^ ^ and they must be 
expressed in terms of the fine-grained vorticity distribution p(r, a). This implies that E and F can be computed at any time from 
the coarse-grained flow while the r^>i are only accessible from the initial conditions or from the fine-grained flow. Since they arc 
not accessible from the macroscopic dynamics, they behave as "hidden constraints". 
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In the case of flows that are forced at smaU-scales, EUis et aJ. [22] and Chavanis [23] have proposed to treat 
the constraints associated with the fragile moments r„>i canonically and fix the Lagrange multipliers a„>i, 
or the prior distribution xC^); instead of r„>i (by contrast, the robust constraints F and E are still treated 
microcanonically) . If wc view the vorticity levels as species of particles, this amounts to fixing the chemical 
potentials instead of the total number of particles in each species. It is argued that the prior x(o") is imposed by 
the small-scale forcing so it must be regarded as given a priori (it is external to the system under consideration). 
In this point of view, one must work with the relative entropy 



pin 



drda, (23) 



where x(o") is assumed given. The relative entropy (23) can be viewed as the Legendre transform Sx = S — 
X)n>i '^nTn of the mixing entropy (19) with respect to the fragile moments r„>i = J uj^dr [23]. The equilibrium 
state is obtained by maximizing (23) while conserving only the robust constraints E and F. This again yields the 
Gibbs state (20) except that now x{c) is fixed a priori by the small-scale forcing. The optimal coarse-grained 
vorticity is given by (21) where -F(<I>) — — (Inx)'('I') is now completely specified by the prior xi'^)- Note that uJ 
satisfies a result of large deviations [22] . It maximizes a generalized entropy [23] : 



sp] = -j 



C(aJ)dr, (24) 



at fixed circulation and energy, where C7 is a convex function directly determined by the prior according to the 
relation [24]: 



UJ UJ 

C{Uj) = - I F-''{x)dx = - I [{\-a.x)']-''{-x)dx. (25) 

In this context, Chavanis [3,23] has proposed a thermodynamical small-scale parameterization of 2D forced geo- 
physical flows where the coarse-grained vorticity evolves according to 



^ + u-VTd=S/ ■ {D 



''--Mi''* 



10 = -Alp, (26) 



r DVlJ ■ Vthdr 1 , , 

J C" \uj) 

This relaxation equation for ZU(r, €) conserves only the robust constraints (circulation and energy) and increases 
the generalized entropy (24) fixed by the prior vorticity distribution xC^)- It differs from the relaxation equations 
(22) of Robert & Sommcria in the sense that the specification of the prior xip) (determined by the small-scale 
forcing) replaces the specification of the Ceisimirs (determined by the initial conditions). However, in both models, 
the robust constraints E and V are treated microcanonically (i.e. they are rigorously conserved). We thus have 
to solve one differential equation (26) instead of N coupled equations (22) for each discretized level a^- The 
relaxation equation (26) enters in the class of generalized mean-field Fokker-Planck equations presented in Sec. 2. 
It is however associated with a microcanonical description since /3(t) varies in time so as to conserve the energy 
£[3]. 

Because of incomplete relaxation, the coarse-grained vorticity can converge towards a metaequilibrium state 
which is different from that predicted by the statistical theory [11]. This coarse-grained vorticity field is a stable 
stationary solution of the 2D Euler equation which is incompletely mixed. The 2D Euler equation admits an 
infinite number of stationary solutions specified by an arbitrary ui — -ip relationship. A solution of the form uJ = 
/(V'), where / is monotonic, which maximizes an i?-function H\U}\ = — J C{ZJ)dr (where C is convex) at fixed 
energy E and circulation F is nonlinearly dynamically stable with respect to the 2D Euler equation [22,23]. For 
example, Tsallis functional Hq = — f (i^' — cJ)dr is a particular //-function. Its maximization at fixed E and F 
determines a particular class of stationary solutions of the 2D Euler equations (polytropic vortices) corresponding 
to aJ = [A — l3{q — l)^p/qY^^''~^\ Minus the enstrophy F2 = / oJ^dr can also be viewed as an //-function leading to 
a linear uj — tp relationship. This is a particular case of the Tsallis functional corresponding to g = 2. Note that the 
relaxation equations (26) and (27) can serve as a numericMl algorithm to construct nonlinearly dynamically stable 
stationary solutions of the 2D Euler equation since the steady state of these equations maximizes, by construction, 
an //-function specified by C(cJ) at fixed circulation and energy. 
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4. Collisionless stellar systems 



For most stellar systems, the encounters between stars are negligible [25] so that the galactic dynamics is 
described by the Vlasov-Poisson system 



df df „^ df 
dt dr dv 



0, 



(28) 



The Vlasov-Poisson and the 2D Euler-Poisson systems present deep similarities as shown in [7]. In particular, 
the Vlasov-Poisson system develops a complicated mixing process in phase space associated with the damped 
oscillations of a proto-galaxy initially out-of-equilibrium and in search of a steady (virialized) state. Due to 
"phase mixing" , a collisionless stellar system can reach an "equilibrium" state (metaequilibrium) on a very short 
timescale. This collisionless relaxation has been called violent relaxation. Lynden-Bell [26] proposed a statistical 
mechanics of this process which is similar to that exposed previously in 2D turbulence. However, this process 
takes place in a six-dimensional phase space instead of the plane [9]. 

In the two-levels approximation where f — rjo or f — 0, the optimal probability p{r, v) of finding the level r]o 
in r, V at equilibrium is obtained by maximizing the mixing entropy 



S\p] = -J {plnp+ (1 -p)ln(l -p)}drdv, 



(29) 



at fixed 7o = / pdrdv (or total mass M) and energy E = ^ J fv^dvdv + | / p$dr. Using the method of Lagrange 
multipliers, the optimal probability p*(r, v), and consequently the coarse-grained distribution function / = p^rfo, 
is given by 



/(r,v) 



1 + Ae''''o(^+*) 



(30) 



Apart from a reinterpretation of the constants, this is similar to the Fermi-Dirac distribution [26]. Using the 
analogy with 2D turbulence, Chavanis et al. [9] have proposed a dynamical model to describe the relax:ation of 
the system towards equilibrium. By using the MEPP, they obtain in the two-levels case an equation of the form 



^ -I- V — - V$ ■ — = — 
dt dv dv dv 



D 



%+mfino 



/)v 



A$ = 47rG 



jldv, 



(31) 



m = - 



/_D|^.vdrdv 
jDf{no-f)v^drdv' 



(32) 



This relaxation equation enters in the class of generalized mean-field Fokker-Planck equations of Sec. 2 associated 
with the Fermi-Dirac entropy in phase space S[f] = — J{{f /vo) ln(//??o) + (1 — f /Vo) lii(l — //'?o)}drdv and with 
a Newtonian potential u = —G/\r — r'j. Note, however, that in the present case the inverse temperature evolves 
with time according to Eq. (32) so as to conserve the energy. Therefore, Eq. (31) increases the entropy S at fixed 
E and M. This describes a microcanonical situation [3]. 

In the general case, we get a system of relaxation equations for p(r, v,7j,t) of the form [9]: 



dp 
dt 



+ v 



dp 
dr 



dp 

dv 



d_ 
dv 



D 



dp 

dv 



+mpiri-f)v 



(33) 



They converge at equilibrium towards the Gibbs state 



P(r,v,?7) 



-(/3e-Hc«)r) 



Z{r,vy 

which maximizes the mixing entropy 



p In p dvdvdrj, 



(34) 



(35) 
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at fixed E, M and 7(7?) = J pdvdv. At statistical equilibrium, the coarse-grained distribution function f = J pridr) 
is given by 

We note that the statistical theory of violent relaxation leads to non-standard (non-Boltzmannian) coarse-grained 
distribution functions. This is due to the existence of an infinite family of conserved quantities, the Casimirs, which 
play the role of "hidden constraints" in our general interpretation of non-standard distributions. The coarse-grained 
distribution functions arising in theories of violent relaxation can be viewed as sorts of superstatistics [24] since 
they are expressed as superpositions of Boltzmann distributions weighted by a non-universal factor depending on 
the initial conditions. 

Like in 2D turbulence, the violent relaxation of coUisionlcss stellar systems is usually incomplete [26]. The 
Vlasov equation admits an infinite number of stationary solutions specified by the Jeans theorem. A solution of 
the form / = /(e) with /'(e) < which maximizes an //"-function H[f] = — J C(/)drdv (where C is convex) at 
fixed energy E and mass M is nonlinearly dynamically stable with respect to the Vlasov-Poisson system [27]. This 
determines a subclass of spherical stellar systems. For example, the Tsallis functional Hq = — ^zrr / (/^ — f)drdv 
is a particular //-function. Its maximization at fixed E and M determines a particular class of stationary solutions 
of the Vlasov-Poisson system with DF / = [A — /3(g — l)e/q\^^^''~^^ called stellar polytropes in astrophysics (the 
standard polytropic index is related to the q parameter by n = 3/2 -|- l/{q — 1)). Like in 2D turbulence, we can 
use generalized Fokker-Planck equations as numerical algorithms to construct stable stationary solutions of the 
Vlasov equation (see [3]). 



5. Self-gravitating Brownian particles and related systems 

In a series of papers, Chavanis & Sire [15] have introduced and studied a model of self-gravitating Brownian 
particles. In this model, the particles interact via self-gravity but, in addition, they experience a friction force 
(against an inert medium) and a stochastic force. This Brownian model is described by the canonical ensemble in 
contrast to the usual Hamiltonian model of stellar systems which is described by the microcanonical ensemble. In 
the mean-field limit, the evolution of the distribution function of the Brownian gas is governed by the Kramers- 
Poisson system. In the strong friction limit, we get the Smoluchowski-Poisson (SP) system 



dp 
dt 



|(TVp-FpV$) 



A$ = SdGp. 



(37) 



(38) 



Comparing with Eq. (6) , this corresponds to an isothermal equation of state p — pT associated with the Boltzmann 
entropy S = — J plnpdr. In d > 2, there exists a critical temperature Tc for box-confined systems. For T > Tc, 
the system reaches a stable equilibrium state corresponding to the Boltzmann statistics 



Peg = Ae 



-/3* 



(39) 



For T < Tc, there is no equilibrium state anymore. The density blows up and finally forms a Dirac peak (conden- 
sate) [15]. Chavanis & Sire [16] have also studied more general models where the diffusion of particles is anomalous 
so that the evolution of the particles is described by the nonlinear Smoluchowski-Poisson (NSP) system 



dt 



-{KVp-' + pW^) 



A$ = SdGp. 



(40) 



(41) 



Comparing with Eq. (6), this corresponds to a polytropic equation of state p = Kp'^ associated with the Tsallis 
entropy S = — ^ / {p'^ — p)dr. The steady states (polytropes) reproduce the statistics introduced by Tsallis 



Peq 



(42) 
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The nonlinear Smoluchowski equation (40) also arises in the physics of porous media where anomalous diffusion is 
due to the complicated (fractal, multi-fractal) phase-space structure of the medium and $ is an external potential 
[5] . More generally, we can consider the evolution of particles described by the generalized Smoluchowski- Poisson 
(GSP) system 



dt 



-(Vp + pV$) 



A$ = SdGp. 



(43) 



(44) 



The case of self-gravitating Brownian fermions where p = p(p) is the Fermi-Dirac equation of state has been 
studied in [17]. Some general properties of the GSP system, including the Virial theorem and a stability analysis 
of the solutions have been obtained in [18]. 

Wo can also consider other forms of potential of interaction than the gravitational one. For example, the 
Smoluchowski-Poisson system with the repulsive Coulombian potential describes electrolytes in the theory of 
Debye-Hiickel [28]. On the other hand, when the particles evolve on a ring and the potential of interaction is 
truncated to one Fourier mode, we get the Brownian Mean Field (BMF) model [19]: 



dt 39^ 2TrddYj 



\{e-e')p{e',i)de' 



(45) 



which is the canonical version of the Hamiltonian Mean Field (HMF) model. 

Finally, a dynamical model of the Bose-Einstein condensation has been studied recently by Sopik et al. [20]. 
The dynamical evolution of a gas of non-interacting bosons in contact with a thermal bath is governed by the 
bosonic Kramers equation 



iv 



TVkP + p(l + p)k 



(46) 



which takes into account an inclusion principle [4] . This equation enters in the generalized class of Fokker-Planck 
equations of Sec. 2. It is associated with the Bose-Einstein entropy S = — /{plnp — (1 -|- p)ln(l -|- p)}dk. We 
have found that the Bose-Einstein condensation (in phase space) shares deep similarities with the collapse of the 
self-gravitating Brownian gas (in position space). In particular, for d > 3, there exists a critical temperature Tc. 
The density blows up for T < Tc leading ultimately to a Bose condensate in k-space [20]. 



6. Chemotactic aggregation of bacterial populations 

The name ehemotaxis refers to the motion of organisms (amoeba) induced by chemical signals (aerasin). In 
some cases, the biological organisms secrete a substance that has an attractive effect on the organisms themselves. 
Therefore, in addition to their diffusive motion, they move systematically along the gradient of concentration of 
the chemical they secrete (chemotactic flux). When attraction prevails over diffusion, the ehemotaxis can trigger 
a self-accelerating process until a point at which aggregation takes place. This is the case for the slime mold 
Dictyostelium Discoideum and for the bacteria Escherichia coli. 

A model of slime mold aggregation has been introduced by Keller & Segel [6] in the form of two PDE's: 

^ = V • (DaVp) - V • (i3iVc), (47) 

dc 

■^ = -k{c)c+fic)p + DcAc. (48) 

In these equations p(r, t) is the concentration of amoebae and c(r, t) is the concentration of aerasin. Aerasin is 
produced by the amoebae at a rate /(c). It can also be degraded at a rate fc(c). Aerasin diffuse according to Pick's 

law with a diffusion eoeffieient D^- Amoebae concentration changes as a result of an oriented chemotactic motion 
in the direction of a positive gradient of aerasin and a random motion analogous to diffusion. In Eq. (47), -D2(p, c) 
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is the diffusion coefficient of tlie amoebae and Di (p, c) is a measure of the strength of the inffuence of the acrasin 
gradient on the flow of amoebae. This chemotactic drift is the fundamental process in the problem. 
A first simplification of the Keller-Segel model is provided by the system of equations 

dp 
at 

dc 



= £)Ap-xV-(pVc), 



(49) 



dt 



= D'Ac + ap — be, 



(50) 



where the parameters are positive constants. An additional simplification, introduced by Jager & Luckhaus [29], 
consists in ignoring the temporal derivative in Eq. (50) . This is valid in the case where the diffusion coefficient D' 
is large. Taking also 6 = 0, we obtain 



^=7?Ap-xV-(pVc), 



Ac: 



-Ap, 



(51) 
(52) 



where A = a/D' . These equations are isomorphic to the Smoluchowski-Poisson system (37)-(38) which describes 
self-gravitating Brownian particles in a high friction limit. The analogy between self-gravitating systems and the 
chemotaxis of bacterial colonies is developed in [17]. 

The KcUcr-Segel model ignores clumping and sticking effects. However, at the late stages of the blow-up, when 
the density of amoebae has reached high values, finite size effects and stickiness must clearly be taken into account. 
As a first step, we have proposed in [17] to replace the classical equation (51) by an equation of the form 



• xV • (p(l - p/fTo)Vc), 



(53) 



which enforces a limitation p < ctq on the maximum concentration of bacteria in physical space. This equation 
decreases the Lyapunov functional 



F[P] 



'_£_ 


ln^ + 






(To 


(To 


0(^J 





dr. 



(54) 



This functional can be interpreted as a free energy F = E — T^ffS associated with a Fermi-Dirac entropy in 
physical space [3,17] where Tg// = D/x is an effective temperature. This form of entropy can be obtained by 
introducing a lattice model preventing two particles to be on the same site. The lattice creates an exclusion 
principe in physical space similar to the Pauli exclusion principle in phase space for fcrmions. Then, S[p\ can bo 
obtained by a standard combinatorial analysis respecting this exclusion principle. The equilibrium states of Eq. 
(53) are given by a Fermi-type distribution in physical space 



P = 



l + Ae^ 



(A > 0), 



(55) 



which minimizes the effective free energy (54) at fixed mass. On the other hand, the relation between the con- 
centration of amoebae and acrasin may be more complex that simply given by the Poisson equation (52). For 
example, taking 6 7^ in the original Keller-Segel model, we obtain a relation of the form 



Ac- 



-Ap, 



(56) 



where k% = b/D' . The second term is similar to the Dcbye shielding in plasma physics or to the Rossby shielding 
in Quasi-Geostrophic (QG) fiows. We note that Eq. (53) enters in the class of nonlinear mean-field Fokker-Planck 
equations of Sec. 2 associated with a Fermi-Dirac entropy in physical space. We also note the resemblance with the 
relEixation equations (17) of 2D turbulence except that here /3 = x/D is fixed (which is similar to a "canonical" 
situation) while in 2D turbulence E is fixed (microcanonical situation). 

More generally, in the primitive Keller-Segel model (47)- (48) the coefficients Di and D2 can both depend on 
the concentration of particles. This takes into account "hidden constraints" that act on the microscopic dynamics. 
Equation (53) is an example of nonlinear mean-field Fokkcr-Planck equations where the mobility depends on the 
concentration while the diffusion coefficient is constant. Alternatively, if the diffusion coefficient depends on the 
concentration while the mobility is constant, we can write Eq. (47) in the form 



f^ = V.[x(Vp-pVc)], 



(57) 
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where p — p{p) is an effective pressure. Tfierefore, the Keller-Segel model is an example of nonlinear mean-field 
Fokker-Planck equations presented in Sec. 2. We can thus develop an Effective Thermodynamical Formalism 
(ETF) to investigate the chemotactic problem. We can also introduce an inertial model of chemotaxis [17,18] 
based on the hydrodynamic equations 

^+V.{pu)=0, (58) 
^ + (u-V)u = -ivp + Vc-Cu, (59) 

p^=-k{c)c+f{c)p + D,Ac. (60) 

These hydrodynamical equations can be derived from generalized Fokker-Planck equations of the form (1) [18]. 
They are similar to the damped Euler equations (9)- (10) for self-gravitating Brownian particles [3] except that 
the Poisson equation is replaced by a more general field equation taking into account the specificities of the 
chemotactic model. For ^ = 0, we recover the hyperbolic model introduced by Gamba et al. [30] and for ^ +oo, 
we can neglect the inertial term in Eq. (59), leading to pu — — |-(Vp + pV$) + 0(^~^). When substituted in Eq. 
(58), we recover the parabolic Keller-Segel model (57). It has been shown that parabolic models of chemotaxis 
lead to point-wise blow-up while hyperbolic models show the formation of a network that is interpreted as the 
beginning of a vasculature. This phenomenon is responsible of angiogenesis, a major factor for the growth of 
tumors. 



7. Conclusion 



In this paper, we have given several physical examples where the evolution of the system is governed by 
nonlinear mean-field Fokker-Planck equations introduced in [3]. These equations are consistent with an Effective 
Thermodynamical Formalism (ETF) and lead to non-standard distributions at equilibrium. We stress that: (i) 
Tsallis distributions and Tsallis generalized entropies just represent a particular case of this general formalism, 
(ii) The physical reason why non-standard distributions arise in a problem is different from case to case and 
must be sought for each particular system. In general, this reflects the existence of "hidden constraints" [3] that 
can be of various forms as discussed in the Introduction. If the system does not "mix well" (for some reason) 
there is no universal form of entropy to account for non-ergodicity. (iii) In many cases, the resemblance with a 
"generalized thermodynamical formalism" is essentially effective or formal. In some cases, there is no relation 
with thermodynamics at all. Yet, it may be of interest to develop a thermodynamical analogy and use the same 
vocabulary as in thermodynamics. This allows to unify the formalisms and transpose the results obtained in one 
context to another context. Indeed, we have found many analogies between systems of a very different nature: 
stellar systems, vortices, bacteria, Bose-Einstein condensation... [3,7,17,20]. 
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